Probing the photonic local density of states with electron energy loss spectroscopy 



F. J. Garcia de Abajo^ and M. Kociak^ 

' Instituto de Optica - CSIC, Serrano 121, 28006 Madrid, Spain 
^ Laboratoire de Physique des Solides, CNRS, Universite Paris Sud XI, F 91405 Orsay, France 

(Dated: February 2, 2008) 

Electron energy-loss spectroscopy (EELS) performed in transmission electron microscopes is 
shown to directly render the photonic local density of states (LDOS) with unprecedented spatial 
resolution, currently below the nanometer. Two special cases are discussed in detail: (i) 2D photonic 
structures with the electrons moving along the translational axis of symmetry and (ii) quasi-planar 
plasmonic structures under normal incidence. Nanophotonics in general and plasmonics in particu- 
lar should benefit from these results connecting the unmatched spatial resolution of EELS with its 
ability to probe basic optical properties like the photonic LDOS. 
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While a plethora of nanophotonic structures are cur- 
rently being devised for diverse applications like achiev- 
ing single molecule sensitivity in biosensing [H or molding 
the flow of light over nanoscale distances for signal pro- 
cessing 0, no optical characterization technique exists 
that can render spectroscopic details with truly nanome- 
ter spatial resolution. The need for that kind of tech- 
nique is particularly acute in nanometric plasmonic de- 
signs that benefit from sharp edges and metallic surfaces 
in close proximity to yield large enhancements of the elec- 
tromagnetic field. 

Scanning transmission electron microscopes (STEM) 
can plausibly cover this gap, as they perform electron 
energy loss spectroscopy (EELS) with increasingly im- 
proved energy resolution that is quickly approaching the 
width of plasmon excitations in noble metals J3| and with 
spatial resolution well below the nanometer [J| . The con- 
nection between EELS and photonics can be readily es- 
tablished when low-energy losses in the sub-eV to a few 
eV range are considered, compatible with typical pho- 
ton energies in photonic devices. A formidable amount 
of information is available in the literature for this so- 
called valence EELS, including for instance studies of 
single nanoparticles of various shapes [1, Q , interacting 
nanoparticles Q , thin films 0] , composite metamaterials 
0, and carbon nanostructures 0. Many of these re- 
ports are relevant to current nanophotonics research, in 
which the optical response of nanoparticles, nanoparticle 
assemblies, and patterned nanostructures plays a central 
role. In this context, EELS has been recently demon- 
strated to image plasmon modes with spatial resolution 
better than a hundredth of the wavelength in triangu- 
lar nanoprisms i^. However, despite significant progress 
from the theoretical side [l^ , no synthetic and universal 
picture has emerged to explain the spatial modulation of 
EELS measurements on arbitrary nanostructures. 

In this Letter, we show that EELS provides direct in- 
formation on the photonic local density of states (LDOS), 
and thus it constitutes a suitable tool for truly nanomet- 
ric characterization of photonic nanostructures. A rigor- 



ous derivation of this statement is offered, illustrated by 
numerical examples for both translationally-invariant ge- 
ometries and planar structures. Our results allow directly 
interpreting EELS data in terms of local photonic proper- 
ties that encompass the full display of optical phenomena 
exhibited by nanostructures, ranging from localized and 
propagating plasmons in patterned metallic surfaces Q 
to band gaps in dielectric photonic crystals 

Green tensor and LDOS.- The optical response of a 
nanostructure is fully captured in its electric Green ten- 



sor and its LDOS [1^. In particular, the electric field 



produced by an external current density j(r, cj) in an in- 
homogeneous medium of permittivity e(r, w) can be writ- 
ten in frequency space lu as 



E(r,cj) = -iniuj J dr'G(r, r', w)j(r', w) 



(1) 



in terms of G, the electric Green tensor of Maxwell's 
equations in Gaussian units, satisfying 



V X V X G(r,r',cj) - {uj^ /c^)e{r,uj)G{r,r' ,uj) 
= 4'5(r-r') 



(2) 



and vanishing far away from the sources. 

We then define the LDOS projected along unit a vector 
n as Il3l 



-2w 



Im{n ■ G(r, r, uj) ■ n}. 



(3) 



In free space, the uniform LDOS is known from black- 
body theory: 



(4) 



Similar to its electron counterpart in solid state physics, 
the photonic LDOS equals the combined local intensity 
of all eigenmodes of the system under consideration, pro- 
vided they are well defined (e.g., in the absence of lossy 
media) |1J|. An alternative interpretation, which holds 
even in the presence of lossy materials [l5l |. comes from 



2 



the realization that (47r^a;D^/?t) p is the decay rate for an 
excitation dipole strength D [lai II3l ■ Finally, we point 
out that a complete definition of the LDOS should in- 
clude a magnetic part [l^, which is however uncoupled 
to our fast electrons. 

Energy loss probability.- The energy loss suffered by 
a fast electron passing near an inhomogeneous sample 
and moving with constant velocity v along a straight line 
trajectory r ~ re(t) can be related to the force exerted by 
the induced electric field E'""^ acting back on the electron 
as [lil 



(5) 



AE = e J dtv- E'"'^[re(i), t] = hcoduj T{, 

where the — e electron charge has been included (i.e., 
AE > 0) and 

r(t^) = / dtRe\e-'"^'w ■-E""^\rJt),uj]] (6) 
nhuj J ^ ' 

is the loss probability. The Fourier transform 

(7) 



duj 



E'"'^(r,i)= / — c-"^*E'"'^(r,t^) 
J 27r 

has been introduced and the property E'"'^(r,CL') = 
[E™'*(r, -oj)]* has been used. 

The external current density corresponding to the 
moving electron is now given by 



j(r,c^) = -ev J dte'^'d[r-r,{t)]. 



(8) 



Assuming without loss of generality that the velocity vec- 
tor is directed along the positive z axis and using the 
notation r = (R, z), with R = (x, y), the current density 
reduces to 



j(r,w) = -e<5(R-Ro)e" 



z/v , 



(9) 



where Rq = {xo,yo) is the 2D impact parameter of the 
electron trajectory relative to the z axis. Inserting Eq. 
© into Eq. ([1]), and this in turn into Eq. ©, we find 



r(Ro,c^) 



(10) 



dtdt 



where Gzz — z ■ G ■ z, the dependence of the loss proba- 
bility r on Rq is shown explicitly, and G""'^ denotes the 
induced Green tensor obtained from G by subtracting the 
free-space Green tensor. 

Relation between EELS and LDOS.- Noticing that 
Ze{t) — vt, the time integrals of Eq. pH)) yield the Fourier 
transform of the induced Green tensor with respect to z 
and z', G"'^(R, R', g, — q', cj), in terms of which F be- 
comes 



F(Ro,c.) = -^lm{GL"/(Ro,Ro,g,-g,t^)}, 



Pi.{Ro,q,uj), 



(11) 



where q = uj/v and we have defined 



-2u; 



Im{n-G(R,R,g, -(7,cj)-n} (12) 



as a generalized density of states that is local in real space 
along the R directions and local in momentum space 
along the remaining z direction, parallel to the electron 
velocity vector. 

The value of q — lu/v reflects conservation of energy 
and momentum in the transfer of excitations of frequency 
uj and momentum q from the electron to the sample. It 
is interesting to note that for an electron moving in an 
infinite vacuum one has 



Pz(R.,g,w) = L^—^0{uj/c-q) 



(13) 



which is always zero for subliminal electrons moving with 
velocity v < c. Here, L is the quantization length 
along z. As expected, the vacuum density of states 
does not contribute to the EELS signal, and this al- 
lows dropping the superscript ind from Eq. (fTTjl . as it 
has been implicitly assumed already when writing Eq. 
([T^. Incidentally, the Cherenkov effect is deduced from 
/9z(R.,g,w) = £(cj/27rc2)(l - k^/q^e)0{ew^/c-q^), valid 
for a homogeneous dielectric of real permittivity e. 

The proposed formalism, yields exactly the same re- 
sults as any other local, retarded theory, but it provides 
a new paradigm for understanding EELS as connected 
to a local quantity: the LDOS. This is in contrast to 
the common view of EELS as a tool capable of retriev- 
ing local electronic properties hidden in the permittivity, 
which is in turn involved in the non-local response of in- 
homogeneous structures probed by the electrons and has 
resulted in endless discussions regarding how to eliminate 
delocalization. 

2D systems possessing translational invariance.- If the 
sample under consideration is translationally-invariant 
along z, the Green tensor G(r,r',ti;) depends on z and z' 
only via z — z', and thus one can write 



G{v,v\uj) = J ^G{R,R',q,co)e 



iq(z-z') 



(14) 



where G(R,R',g,w) = (1/L)G(R, R', g, -g, w). Accord- 
ingly, the local density of states can be decomposed into 
momenta components q along the z axis, 



(15) 



where 



Pfi(R, g,cj) = — ^Im{n • G(R, R, g,w) • n} 

TT 

= j-Ph(R,q,uj)- 



(16) 



Finally, combining Eqs. pTjl and one finds a relation 
between the loss probability per unit of path length and 



the LDOS, 



Pi{Ro,q,uj), 



with q = Lo/v. 

Eq. (|17p provides a solid hnk between LDOf 
EELS in a wide class of geometries that include 
trajectories in semi-infinite surfaces and thin films 
terestingly, only q> uj/c values are probed, lying o\ 
the light cone, and therefore difficult to study via 
cal techniques. Trapped modes such as surface-pla 
polaritons lie in that region and are a natural targ 
application of our results. Besides, the present stud 
be directly applied to cathodoluminescence (CL) i 
dielectric structures, in which energy loss and CL 
sion probabilities are identical. 

A connection between the photonic density of f 
in the momentum space and EELS has been previ 
reported for electrons moving parallel to pores in 21 
assembled alumina photonic crystals pjj. Howeve 
above derivation is the first prove to our knowledge 
a formal relation exists between LDOS and EELS. 

For illustration, we offer in Fig.[l]the EELS proba 
for electrons moving inside a finite hexagonal 2D c: 
of aligned Si nanowires, calculated with the bou] 
element method (BEM) [23|. The photon wave! 
range under consideration includes two Mie modes i 
isolated wire (broken curve), and significant hybr 
tion between neighboring wires takes place in the ; 
The loss probability varies relatively smoothly wit 
pact parameter, a behavior which was expected i 
LDOS for photon wavelengths relatively large com 
to the cylinders diameter. Interestingly, the loss p 
bility takes significant values in the intersticial regions, 
several tens of nanometers away from the Si. Finite struc- 
tures like that considered in Fig. [T] exhibit a colorful evo- 
lution of modes, the analysis of which can be useful for 
instance in the design of microlaser cavities. In infinite 
crystals, the loss probability exhibits Van Hove singular- 
ities 2l|, which follow rigorously those of the LDOS in 
translationally-invariant systems according to Eq. (fT7|) . 

Planar geometries.- As microchip features continue to 
shrink, lithographically-patterned metal structures are 
becoming natural candidates to replace current electronic 
microcircuits. The new structures will operate at fre- 
quencies above the THz, rather than GHz, and will carry 
electric signals strongly mixed with the electromagnetic 
fields that they generate in what is known as surface plas- 
mons. The optical properties of metallic planar struc- 
tures are routinely obtained using scanning near-field op- 
tical microscopy (SNOM), although the lateral resolution 
of this technique can hardly reach 50 nm. Here again 
EELS renders much higher lateral resolution (down to 
the nm) and permits obtaining information directly re- 
lated to the LDOS. Besides the formal relation between 
EELS and the momentum-resolved LDOS expressed in 
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FIG. 1: EELS probability for 200-keV electrons moving par- 
allel to an array of 13 Si nanowires of 300 nm in diameter 
(solid curve), as compared to an isolated nanowire (broken 
curve). The nanowires are arranged in a hexagonal lattice 
with a period of 350 nm. The contour-plot insets show the 
impact-parameter dependence of the EELS probability for se- 
lected energy losses (labeled A-G), with the position of the 
electron beam considered in the curves indicated by an open 
symbol and the cylinder contours shown by dashed curves. 



Eq. (|12|). we offer in Fig. [2] a more detailed comparison 
between EELS and p(r, uj) (local in all spatial directions) 
for an Ag disk calculated with BEM |20|. This figure 
proves how the EELS probability can mimic quite closely 
Pzi^, w), with z perpendicular to the planar structure and 
parallel to the electron velocity. This resemblance holds 
for different accelerating voltages, making the interpre- 
tation robust with respect to experimental details. Fig. 
[5] provides a solid example supporting the use of EELS 
to measure plasmon intensities with unprecedented lat- 
eral resolution. A first demonstration of such type of 
measurements has been recently reported [3]. 

We thus conclude that the energy loss probability is 
directly related to the local density of states in arbitrary 
systems, where we understand locality in real space for 
the directions perpendicular to the electron trajectory 
and in momentum space along the direction of the elec- 
tron velocity vector. In 2D systems and for electrons 
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FIG. 2: Relation between EELS and LDOS in planar geometries, (a) We consider an Ag disk of height 10 nm and radius 30 nm. 
The electrons move along z, perpendicular to the disk. The EELS probability for 200-keV (d) and 100-keV (e) electrons mimics 
closely the z-projected LDOS in a plane 10 nm above the disk (b), and less closely the unprojected LDOS {p = Px + Py + pz) 
(c). The cathodoluminescence emission only picks up part of the inelastic signal (f). 



moving along the direction of translational symmetry, the 
loss probability is exactly proportional to the photonic 
local density of states projected on the trajectory and 
decomposed into parallel momentum transfers q. Nu- 
merical examples have been presented showing a similar 
relation between LDOS outside planar metallic disks and 
EELS spectra for electrons traversing them perpendicu- 
larly. Our results provide a solid foundation for the use 
of EELS performed in STEMs to directly probe photonic 
properties of nanostructures. 
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